On the Salecker-Wigner-Peres clock and double barrier tunneling 
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In this work we revisit the Salecker-Wigner-Peres clock formalism and show that it can be directly 
applied to the phenomenon of tunneling. Then we apply this formalism to the determination of the 
tunneling time of a non relativistic wavepacket, sharply concentrated around a tunneling energy, 
incident on a symmetric double barrier potential. In order to deepen the discussion about the 
generalized Hartmann effect, we consider the case in which the clock runs only when the particle 
can be found inside the region between the barriers and show that, whenever the probability to find 
the particle in this region is non negligible, the corresponding time (which in this case turns out to 
be a dwell time) increases with the barrier spacing. 

PACS numbers: 03.65.Xp, 73.40.Gk 
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I. INTRODUCTION 

An unambiguous definition of a tunneling time is an 
important problem in quantum mechanics, due to both 
its applications and its relevance to the foundations of 
the theory. It is, however, a problem which has eluded 
physicists since the beginnings of the quantum theory. 
Many attempts have been made to define such a time 
scale (see [H, 0, [3, 0] for reviews). However, most of these 
definitions (phase time, dwell time, Larmor time, etc.), 
while being valid to describe some specific characteristic 
of the tunneling process, also present difficulties if one 
tries to interpret them as traversal times in general. 

Perhaps the most striking of the above mentioned diffi- 
culties concerns the issue of superluminality, a direct con- 
sequence of the Hartmann effect Q, which asserts that 
the phase time saturates for opaque barriers. More re- 
cently, the Hartmann effect has been considered for the 
double barrier potential, and it was verified that in the 
opaque limit the phase time does not depend on the 
spacing between the barriers either, a phenomenon re- 
ferred to as the generalized Hartman effect |6[ (also see 
0,113). Altho ugh there are no real paradoxes associ- 
ated with these phenomena, since it is well known that 
the group velocity cannot be associated to the signal ve- 
locity in such a situation [lol ]. the subject has originated 
intense debate in the literature (see, for example, [ll| and 
the references there cited). Independently of these con- 
troversial interpretations, it is in fact a counterintuitive 
result that in the opaque limit the phase time does not 
depend on the spacing between the barriers, because one 
could expect group velocity to have the usual physical 
meaning in that free region. Subsequent investigations, 
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both in the non-relativistic [l2| and in the relativistic [T3[ 
cases, indicate that this lack of dependence may be an ar- 
tifact of the opaque limit, but the subject still deserves 
further investigation. 

A fruitful avenue of investigation on tunneling times 
considers the use of quantum clocks. A quantum clock 
is a secondary dynamical system weakly coupled to the 
system of interest and having a degree of freedom evolv- 
ing uniformly in time. One of the most prominent works 
along this line leads to the Larmor time fbil Il5j ]. but 
other clocks are possible (see for example, [2( and refer- 
ences there cited) . Here we are particularly interested in 
the clock formalism introduced by Salecker and Wigner 
[IH and later revisited by Peres [13], who used it to in- 
vestigate, among other problems, the time-of-ffight for a 
non-relativistic particle (see also [Hj]). The extension for 
a relativistic particle was later done by Davies . 

In [13] Peres also introduced a "time operator" (not 
canonically conjugate to the clock's Hamiltonian) , whose 
expectation values do not lead to sensible results for the 
tunneling time in the presence of a localized potential, as 
was later shown by Leavens et al [2(J, Hl| ■ To overcome 
such a difficulty, these authors proposed a modification of 
the original Salecker-Wigner-Peres (SWP) formalism by 
the introduction of a calibration procedure. However, in 
his treatment of the time-of-flight problem Peres [13] did 
not use directly such an operator, but defined the time 
given by the clock (t c ) as the derivative of the phase shift 
of the wavefunction with respect to the perturbation po- 
tential. In this work we demonstrate that Peres' original 
approach, contrary to what is usually stated, can be di- 
rectly applied to the tunneling time problem, without 
the need for calibration. In section [TT] we present a brief 
review of the SWP formalism and clarify some impor- 
tant issues related to it. In particular, we present a sim- 
ple proof of the general result that t c (averaged over the 
scattering channels) is exactly equal to the well known 
dwell time (a result obtained by Leavens through the ex- 
pectation values of the Peres' "time operator" only after 
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calibration). In section Hill we apply the SWP formalism 
to the tunneling through a symmetric double barrier and 
analyze the dependence of t c with the spacing between 
the barriers. It must be noticed that although in this 
case the (transmitted or reflected) time resulting from 
the SWP clock is exactly the dwell time, this formalism 
proves to be operationally better suited to address the 
question of independence or not of the tunneling time 
with respect to the barrier spacing in the limit of opaque 
barriers, providing a simpler procedure for the direct cal- 
culation of the time spent by the wavepacket only between 
the two barriers. Therefore, such approach allows us to 
deepen the discussions about the generalized Hartmann 
effect. In Section |IV] we discuss the results and their 
interpretation. In the Appendix we list the explicit ex- 
pressions for some terms appearing in the expressions for 
the times obtained in section HTTl 



II. 



THE SWP CLOCK AND THE TUNNELING 
TIME PROBLEM 



The free SWP clock consists of a quantum rotor, which 
for a Hilbert space of dimension N has Hamiltonian given 

by E3 



H c = ujJ . 



(1) 



with J = —ihJL, to = and r is the clock's resolution. 
The energy eigenstates are u m = ^g^i ^ e [0, 27r), with 
eigenvalues V m — mftio (m = —j, Another conve- 

nient orthonormal basis for the clock's Hilbert space is 

m 



i 



sin 



uikr) 



9 e [0, 2vr) (2) 



Vk{6) V^N sin [§ (9 - ukr)] 

with k = 0, .... N—l. These are the eigenstates of the 
Hermitian operator 



JV-l 



T = krVk , with V k vi = S k ivi, 



(3) 



k=0 



with eigenvalues kr. The above operator plays the role 
of a "time" operator, despite not being canonically con- 
jugate to the Hamiltonian. The motivation for this iden- 
tification is that, for large TV, the wavefunctions v k are 
sharply peaked around 9 = ujkr, with a width and 
their time evolution is given by 



<<xp [ -^Hct 



v k {9) =v k {9-Lut). 



(4) 



Thus, Vk{6) evolves rigidly and uniformly within the in- 
terval [0, 27r). In particular, for large N the peaks trans- 
late from tokr to Lo(kr + t) (mod 2tt). 

It must be noticed that it is only for times t = nr, with 
n an integer, that 

V k (0 - UJUt) = V {k+n) (mod AT) W) 



such that the whole set of eigenfunctions Vk{9) can be 
obtained from any of them (say Vq) through a sequence 
of (discrete) time translations. This fact is at the ori- 
gin of the discrepancies found by Leavens [2(| between 
the intrinsic time t in ([4J and the expectation value of 
(|3|) whenever t is not an integer multiple of t. To over- 
come such a difficulty, Leavens introduced a calibration 
procedure, later revised by Leavens and McKinnon [2l| . 
designed in such a way that after calibration the clock 
times (given as an average over an ensemble of freely 
running clocks) coincided with the intrinsic ones. 

On the other hand, the above properties of the wave- 
functions Vk under time evolution allowed Peres to con- 
sider them as the proper clock's hand, with the clock's 
"reading" given by the angle cut c (the translation of the 
wavefunction's peak). In applying this approach to the 
one-dimensional scattering of a particle of mass /i by 
a localized potential V(z) confined within the region 
< z < L, the clock-system coupling can be designed 
to measure the time the wavepacket spends within an 
arbitrary region Z\ < z < Z2- In this case the Hamiltonian 
for the coupled system is given by [13 



H = H S + P(z)H c , 



(5) 



where H s = 2- + V(z), H c is the clock Hamiltonian 
(fTj) and P{z) is a projection operator into the interval 
(zi, Z2). Let us consider a particle in a stationary state of 
energy E incident from the left (the results remain valid 
for a wavepacket strongly concentrated around E). For 
the initial (free) clock state we choose, following Peres, 
vq(6). Then, assuming that the highest eigenvalue of the 
clock, V m , is negligible when compared to all the relevant 
energy scales in the problem, the final (asymptotic) state 
of the whole system is given by (171 120 ] 



T e lkz va(9-wt^), z> max {z 2 , L} 
R cr tkz v Q (9 - uj tf) , z<min{zi,0} : 



where T and R are the transmission and reflection coeffi- 
cients (which depend on the energy E) for the system in 
the absence of the clock. From the above expressions and 
Q, one identifies the Peres' transmission and reflection 
times, respectively, by 



t T c {E)- 



dm) 



and t?(E)=-h- 



where 



dm) 



4 mh 



n =0 

(6) 

is the phase delay of transmis- 
sion (reflection) in the presence of the clock, and the 
superscript (m) indicates the m-th clock's eigenstate. 
In deriving the above result it is assumed that in the 
vanishingly weak coupling limit 7 1 ( m ) ps \T \ e^T and 

w |i?|e^H m ' [l7|. 
A relation between the dwell time and ^ can be ob- 
tained by following steps similar to those Winful used 
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to derive a relation between the dwell time and the 
phase times 22 j (see also [23]). In order to do this, 
one must realize that for a stationary incident parti- 
cle and when the clock is in its m-th stationary state 
the problem is reduced to the solution of the time- 
independent Schrodinger equation with a localized po- 
tential V (m )(z) = V(z) + V m V(z) [Hill, whose solution 
outside the potential region is given by 

r T (m) e ikz t z > max { Z2 .£} 

^ m \z) = { (7) 
[ e lkz + R ( - m h~ lkz 1 z<min{zi,0} , 

where k = j-^/2/iE. Considering the Schrodinger equa- 
tion with the potential V^ m '(z) and its complex conju- 
gate we obtain, after taking the vanishingly weak cou- 
pling limit V m ^0, 



2/i dz 



dz 

m) 



v m =o 



\ dV m di 



(8) 



where tp denotes the wavefunction after the limit V m — >0. 
Integrating the above expression over the region (zi,z 2 ) 
we obtain ^ f™dz\tp\ 2 = [...] Z2 - [...] Z1 , where [...] cor- 
responds to the term into brackets in the above expres- 
sion, which is constant for all z > Z2 and for all z < z\, 
because V(z) = in those regions. Taking advantage of 
this fact, we can use any value of z > Z2 to compute the 
bracket [...] Z2 . F° r convenience we choose z into the re- 
gion corresponding to the transmitted wave. In the same 
way, we can choose z into the incident/reflection region 
to compute the bracket This procedure, together 

with ©, yields 



ih 2 k 



dz\tp\ 2 = T 



V dV„ 



R 



v m =o 



which can be rewritten as 



JL 

hk 



\T\ 



n j.( m ) 

h d4> T 



4 m) 



Identifying the incident flux j in — the l.h.s. of the 
above expression is the well known expression for the 
dwell time t D 0,0. Finally, from © we obtain 



<D 



\T\H 



2.T 



\R\ 2 t« 



Although the above relation was also obtained by 
Sokolovski et al. [2J] through the use of Feynman's 
path integrals, our proof is worth mentioning due to its 
simplicity. Analogous results were also obtained in the 
framework of weak measurements or through the Lar- 
mor clock formalism, but involving, in general, complex 
times [25|, [2^, [2?J • Relation @, however, involves only 
real times. Leavens and McKinnon [2l[ showed that us- 
ing the "time operator" , eq. Q , such relation can only 
be obtained after applying their calibration procedure. 
As another important point concerning the above rela- 
tion we emphasize that, as long as the transmission and 
reflection times arc defined by ([6]), no interference term 
enter it (the analogous relation involving phase times nec- 
essarily requires such a term [22j). The validity of such 
a relation has been much debated in the literature (see, 
for instance, [l], 0, [28| for different points of view) and we 
hope that the above derivation helps to clarify the fact 
that it all depends on how the transmission and reflection 
times are defined (see also [!, [25| for further discussions) . 

Finally, we note that when the whole potential V(z) + 
VrnP(z) is symmetric the reflected and the transmitted 
phases differ only by a constant [l5| , which leads to tjj = 

= tf [1, [H| . In such a case it must also be noted that 
any of the expressions in ^ constitute an operationally 
simpler way to calculate the exact dwell time. We shall 
take advantage of this fact in the next section. 



III. THE DOUBLE BARRIER TUNNELING 

Let us now consider a particle having a given energy 
E (or a wavepacket sharply concentrated around this en- 
ergy) incident from the left on a symmetric double barrier 
potential, given by 



v m =o V(z) = V {6(z)9(a - z) + Q(z - d - a)6(d + 2a - z)} . 

(10) 

We will consider only the case E < Vq, characterizing 
a tunneling process. Even though we are chiefly con- 
cerned with the case in which the clock runs only if the 
particle is inside the region separating the two barriers 
(a < z < a + d), it is instructive to first consider the 
clock running when the particle is anywhere within the 
potential region (0, 2a+d). The solution of the time- 
independent Schrodinger equation outside the potential 
region in this case is of the form ijT]). with the transmis- 
(9) sion amplitude given by 



= %ik Vm q 2 m e-^ 2a ^ k {2 sin(p m d)(fc 2 + q 2 m ){p 2 m + q 2 n ) - [(k - iq m ) 2 (2p m q m cos(p m d) 

+ Om - <L) MPmd)) e~ 2q ™ a + {q m - ik) 2 (2p m q m cos(p m d) - (p 2 m - q 2 m ) sin(p m d)) e 2 ^ a ] } _1 , (11) 
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where q r , 



v/2/i(Vo + V m - Ej/h and p n 



^/2[i(E — V m )/h. The corresponding phase is 



4 m) 



~(2a + d)k - tan 1 



0n 



with a rn and /3 m defined as 

a m = 2kq m [2p m q m cos(p m d) cosh(2q m a) 
+ ioli ~ V m ) sin(p m d) sinh(2g m a)] ; 



(12) 



(13) 



0m = ~(k 2 + q 2 l )(p 2 n + qf n )sm(p m d) 

{q 2 m - k 2 ) cos(p m d) sinh(2g m a) 
+ (qfn - Pm){lm - k 2 ) sin(p m d) cosh(2g m a). 

From the symmetry of the total potential (including the 
term due to the clock) it follows that tr> = tf = t£ . So, 
using ([6]) we obtain 



where 



#bar 



/' 



hq a 2 + 0£ 



(17) 



is just the (transmitted) Peres' time obtained by allowing 
the clock to run only when the particle passes within any 
of the two barriers (which corresponds to take p m — > 
k in all the above perturbed expressions, before taking 
the derivative in ©). From the proof presented in the 
previous section, together with the symmetry of the total 
potential, it follows that (|15p and (|T7|) are the dwell times 
spent in the regions between and within the two barriers, 
respectively. 

In order to discuss the generalized Hartman effect we 
specialize our results to the opaque limit qa — > oo, in 
which the Peres' time (fT4"|) for the whole potential region 
saturates to the value 



to=ti=- 



h (a 2 



01) 



hi 
k 



hi 
<i 



(14) 



where ao and 0o are obtained from (|13|) by taking the 
limit V m — > (which corresponds to p m — > k and q m — > 
q = i -\/2/i(Vb — E)) . The explicit expressions for ao and 
0o as well as for h\ and hi are given in the Appendix. 

Now, let us consider the case in which the SWP clock 
runs only when the particle can be found inside the region 
between the two barriers, namely in the interval (a, a+d). 
The (transmitted) Peres' time t^ et in this case can be 
obtained from the above results simply by taking the 
limit q m — > q in (p~2|) before taking the derivative in ([6]) , 
which gives 



^bet _ _ ^ 



hi 



hk 



01 



Using the above expression, we can rewrite (| 14[) as 



t D = t h c Ct + t h c a 



(15) 



(16) 



to = ti 



fi 2kq 

qa^oo hq 2 (k 2 + q 2 ) 



(18) 



Then, in the opaque limit the behavior of the Peres' time 
(dwell time) for the whole potential region is analogous to 
that of the phase time, in which it is independent both 
of the barriers width a and the barrier spacing a, and 
we obtain a version of the generalized Hartman effect. 
The above saturated result could also be obtained from 
the non relativistic limit for the dwell time obtained in 
[l3| (in that reference the dwell time was obtained from 
a relation among the pha se and the dwell times involv- 
ing interference terms [231]). However, methods based on 
the phase time, which is an asymptotically extrapolated 
quantity, are not suitable to study the behavior of the 
time the particle spends only inside the region between 
the barriers. 



In the opaque limit qa —> oo, the expression (|15p im- 
mediately yields ij? ct — > 0. A more careful analysis taking 
into account the leading terms in the asymptotic situa- 
tion in which qa is large, but finite, shows that 



ibct 



A^q 2 



-2qa 



{2kd(k 2 



Akq sin 2 (kd) + (k 2 - q 2 ) sin(2fcd)} 



<ja>l h (k 2 + q 2 ) 



{{k 2 - q 2 ) sin(fcd) - 2kqcos(kd)Y 



(19) 



r 



Apart from terms coming from the multiple reflections 
and/or interference at the barriers, this expression clearly 
displays an increasing (almost linear) dependence of t^ ct 
with respect to the barrier spacing d , for any finite qa. 
Fig. [1] shows the behaviors of the three times entering 
expression (|16p . with increasing d. Fig. [TJi concerns a 
large but finite barrier width a. We can observe the 
almost linear increasing of t hct with increasing d. In 



Fig. [TJ) the barrier width is increased three times and 
we can already observe the tendency to saturation (ex- 
cept by the peaks) of the three times: to and i^ ar tends 
to the saturated value (fl8|) . while t^ ot tends to saturate 
to zero. Increasing even more the barrier width makes 
the peaks vanish. It can also be shown that if the clock 
runs only within the first barrier, the corresponding sat- 
urated dwell time in the opaque limit is exactly the same 
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FIG. 1: (Color online) Peres' clock time dependence with 
respect to the distance d between the two barriers. In both 
figures the three clock times appearing in relation (|16p are 
shown, corresponding to the whole potential region (0, 2a + d), 
to the region between the barriers (a, a + d), and to the region 
within the two barriers (0, a) U (a + d, 2a + d). Natural units 
were used (ft = c=l),so that mass and energies are expressed 
in units of the particle mass /i, and distances and times are 
expressed in units of fi~ . In both figures E = 0.01 and 
V = 0.018. (a) a = 10; (b) a = 30. 

as (|18l) . Summarizing, in the opaque limit qa — > oo the 
transmission amplitude beyond the first barrier goes to 
zero, and its associated phase becomes meaningless. This 
fact corroborates the view that the generalized Hartman 
effect, which asserts the independence of the tunneling 
time on d, i s indeed an artifact of the opaque limit (see 
also (Hill). 

IV. CONCLUDING REMARKS 

In this work we revisited the SWP clock formalism, 
which (despite of being based on a thought experiment) is 
useful to understand the fundamental concepts involved 
in the definition of a time scale for tunneling processes. 
We showed that in Peres' version it can be directly ap- 
plied to the problem of tunneling times. We demon- 
strated that departures from the Peres' approach are not 
necessary, even in the case of localized potentials, if one 
focuses, as Peres did, on the time evolution of the eigen- 
functions of his "time operator" , instead of focusing on 
its expectation values, as did Leavens et al plj . 

Using the Peres' approach, and through a simple ex- 
tension of a proof originally designed by Winful in the 
context of phase times [22|, [23[ , we have shown that the 
Peres' times (|6]), when weighted by the transmission and 
reflection probabilities, averages exactly to the dwell time 
(incidentally, for symmetric potentials any of the two ex- 
pressions (JBJ) provides an operationally simpler way to 



calculate the dwell time) . On the other hand, we did not 
address questions regarding the SWP clock's resolution 
in the presence of a localized potential (see [13, HO]) be- 
cause this issue can be addressed in the framework of the 
weak measurement theory , as suggested in [3^, [H| ■ 
Besides that, in this paper each of the time readings as- 
sociated to the clock is equivalent to a dwell time, which 
is a well established time scale 0, [H, H, Q ■ 

We then applied the SWP formalism to the symmet- 
ric double barrier potential, aiming to analyze the so- 
called generalized Hartman effect. We calculated explic- 
itly the dwell time and verified that in the opaque limit 
it does not depend on the barrier separation, confirming 
the emergence of the generalized Hartmann effect also 
in this case (see also [12|])- However, we added a new 
insight into this debated question by allowing the clock 
to run only inside the region between the barriers, and 
taking into consideration the leading terms when the bar- 
rier width is large, but finite, we unambiguously showed 
that the dwell time increases "almost linearly" with the 
barrier spacing d (apart from terms arising from the mul- 
tiple reflections/interference inside this region). The fact 
that such a behavior is modulated by an exponential de- 
cay exp(— 2qa) can be understood by noticing that the 
dwell time is an average over the probability of finding 
the particle in the interest region (and since this probabil- 
ity decays exponentially with qa, so will tp). Therefore, 
whenever the probability to find the particle in the region 
between the barriers is non negligible, the corresponding 
dwell time depends on the barrier spacing. 

All the above considerations reinforce the conclusions 
arrived in [l2j]> in the context of a Fabry-Perot cavity, 
and in [TI|, for the relativistic tunneling through dou- 
ble barriers, that the generalized Hartman effect is just 
a mathematical artifact of the opaque limit: although 
in that limit the transmission phase is well defined and 
finite, it is meaningless since T itself goes to zero. There- 
fore, any time scale defined in terms of such phase (such 
as phase times and, as seen in section II, dwell times) 
also becomes meaningless in this limit: it corresponds to 
the trivial fact that the particle does not penetrate past 
the first barrier, and therefore it makes no sense to as- 
sociate any time duration to its passage (or dwelling) in 
the region between the two barriers. 

APPENDIX 

In this Appendix we present the explicit expressions 
for oto, /3o j hi, an d ^2, which appear in the expressions 
for the Peres' clock times {H}, CES), and (fTT)) : 

«o = 2kq[2kq cos(kd) cosh(2qa) 

+{q 2 - k 2 ) sin(fcd) smh(2qa)] , 

P = -(k 2 +q 2 ) 2 sm(kd) + 2kq(q 2 ~k 2 )cos(kd)smh(2qa) 
+ (q 2 - k 2 ) 2 sin(fcd) cosh(2ga), 
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while hi = ao7i — P0J2 and hi = ao73 — /?o74, with 
I 



71 = -2k(q 2 + k 2 ) sm(kd) - d{q 2 + k 2 ) 2 cos(fcd) + 2q(q 2 - k 2 ) cos(kd) sinh(2ga) 

- 2qkd(q 2 - k 2 ) sin{kd) sinh(2ga) - 2k(q 2 - k 2 ) sm(kd) cosh(2<?a) + d{q 2 - k 2 ) 2 cos{kd) cosh(2<?a); (A.l) 

72 = 2kq {2q cos(kd) cosh(2qa) ~ 2qkd sm(kd) cosh(2qa) 

- 2ksm(kd)sinh{2qa) + d(q 2 - k 2 ) cos(kd) smh(2qa)} (A.2) 

73 = -4q(q 2 + k 2 ) sin(fcd) + 2k('iq 2 - k 2 ) cos(kd) sinh(2ga) + Akqa(q 2 - k 2 ) cos(kd) cosh(2ga); 

+ Aq(q 2 - k 2 ) sin(kd) cosh(2ga) + 2a(q 2 - k 2 ) 2 sm{kd) sinh(2ga) (A.3) 

74 = 2k {4kq cos(kd) cosh(2qa) + (3q 2 — k 2 ) sin(fcc?) sinh(2ga) 

+ Akq 2 a cos(kd) sinh(2ga) + 2qa(q 2 - k 2 ) sin(fcd) cosh(2ga) } . (A. 4) 

I 
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